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Abstract. Mid- infrared images of the Seyfert 1 galaxy 
Mrk 279 obtained with the ISO satellite are presented to- 
gether with the results of a one-year monitoring campaign 
of the 2.5-11.7/tm spectrum. Contemporaneous optical 
photometric and spectrophotometric observations are also 
presented. The galaxy appears as a point-like source at the 
resolution of the ISOCAM instrument (4-5"). The 2.5- 
11. 7 /mi average spectrum of the nucleus in Mrk 279 shows 
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a strong power law continuum with a = —0.80 ± 0.05 
(F„ oc v a ) and weak emission PAH features. The Mrk 279 
spectral energy distribution shows a mid-IR bump, which 
extends from 2 to 15-20 /an. The mid-IR bump is consis- 
tent with thermal emission from dust grains at a distance 
of ^ lOOlt-d. No significant variations of the mid-IR flux 
have been detected during our observing campaign, con- 
sistent with the relatively low amplitude (~ 10 % rms) 
of the optical variability during the campaign. The time 
delay for H/3 line emission in response to the optical con- 
tinuum variations is r ~ 16.7j^'g days, consistent with 
previous measurements. 
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1. Introduction 

According to the unified model of active galactic nuclei 
(AGN), the central massive black hole, its surrounding 
accreting material and the broaddine region (BLR) are 
all embedded within a dusty region, probably a thick 
molecular torus. Along some directions, the dust extinc- 
tion is sufficient to block all UV, optical, and near-IR ra- 
diation originating in the inner components. A review of 
the arguments which led to this picture is presented in 
Wills (1999). The presence of a universal inflection point 
near 1.2 fim in the spectral energy distribution of radio- 
quiet AGN's strongly suggests that the bulk of the near 
IR flux arises from dust thermal emission (e.g., Barvai- 
nis 1987; Sanders et al. 1989). The corresponding color 
temperature, ~ 1,500 K, matches closely the sublimation 
temperature of graphite, the most resilient of the grain 
constituents. The near IR emission can be variable, and 
therefore originates, at least in part, in a compact region. 
Furthermore, in three AGN, the near IR variations have 
been shown to be delayed with respect to the UV-optical 
variations. The time delay corresponds closely to the light- 
travel time to the dust sublimation radius r- m ] measured 
values of r- m are 400 light days for Fairall 9 (Clavel et al. 
1989), 50 light days for NGC 1566 (Baribaud et al. 1992), 
80 light days for NGC 3783 (Glass 1992) and 32 light days 
for Mrk 744 (Nelson, 1996). 

The emerging picture is one where the near to mid-IR 
emission arises from thermal re-radiation of UV and opti- 
cal photons absorbed by the circumnuclear dust. Various 
models for the geometry and location of this dust have 
been proposed, but the exact configuration of the models 
remains unconstrained due to a lack of suitable observa- 
tional data. 

One can use variability as a tool to probe the inner few 
light years of the dusty regions. Reverberation- mapping 
techniques (Blandford & McKee 1982) have been used ex- 
tensively to map the BLR in several AGN, on scales of 
light days to light months, notably by the International 
AGN WatcbQ consortium (Alloin et al. 1994). A similar 
approach can be used to probe the IR-emitting region, i.e. 
the warm dust component within the obscuring material. 
Given UV flux variations of sufficient amplitudes, a mid- 
IR monitoring campaign of sufficiently long duration and 
adequate sampling rate, it may in principle be possible to 
recover the transfer function of the dust. 

The Infrared Space Observatory (ISO; Kessler et al. 
1996) offerred a unique opportunity to carry out such 

For a complete panorama of the AGN Watch data sets, 
result s, and related studies, see the AGN Watch web p age at 
URL jittp: / / www. astronomy.ohio-state.edu/^agn watch / . 



a spectrophotometric monitoring program in the mid- 
IR. The Seyfert 1 galaxy Mrk 279 (z = 0.0294) was 
selected because its celestial position allows an uninter- 
rupted 12-month visibility window for ISO and it has a 
well-documented variability history in the optical (Oster- 
brock & Shuder 1982; Peterson et al. 1985; Maoz et al. 
1990; Stirpe et al. 1994, the UV (Chapman et al. 1985), 
and X-rays (Reichert et al. 1985). Balmer-line time-delays 
(Maoz et al. 1990; Stirpe & de Bruyn 1991; Stirpe et al. 
1994) suggest a BLR size in the range 6 to 12 light days. 
A search for day-to-day variability across the Balmer-line 
profile was unsuccessful (Eracleous & Halpern 1993). No 
far-IR flux variations were detected with IRAS (Edelson 
k Malkan 1987). 

2. IR observations and data reduction 

Mrk 279 was observed with two of the instruments on 
board the ISO satellite: nine narrow-band filter images 
were obtained with ISOCAM (Cesarsky et al. 1996), while 
spectra were recorded with the PHT-S spectrometer, a 
subsystem of the ISOPHOT instrument (Lemke et al. 
1996). The PHT-S spectrometer covers the 2.5-12 /^m 
spectral range at a mean resolution of ^3150kms _1 , 
with a gap between 4.9 and 5.9 /im. Its entrance aperture 
projects to 24" x 24" on the sky. All PHT-S observations 
were carried-out in an identical fashion. The integrating 
amplifiers were reset every 32 s and on-source measure- 
ments were interleaved with background measurements by 
"chopping" on the sky at a frequency of 1/256 Hz. The 
chopper throw was set to 300". For each of the observa- 
tions, total on-source integration time was 2048 s and total 
observing time (including background measurements and 
instrument overheads) was 4236 s. 

The CAM observations were performed in staring 
mode at a magnification of 3 " per pixel and with a gain 
of 2. Nine different filters were used. The unit integration 
time was 2.1s per readout and there were between 72 and 
197 readouts per image, depending on the filter used. The 
particular sequence of filters was chosen to go from high 
to low illumination so as to minimize the detector stabi- 
lization time. Good stabilization was also guaranteed by 
the relatively large number of readouts per exposure. 

The PHT-S observations were made at 16 different 
epochs, from 1996, February 5 to 1997, February 13 (Ta- 
ble 1). The ISOCAM observations of Mrk 279 were all 
carried out on 1996 February 5, contemporaneously with 
the first of the PHT-S spectra. 

Standard procedures from the CAM Interactive Anal- 
ysis (CIA) software package ^ were used for the reduction 
of the ISOCAM data (Ott et al. 1997). The full width at 
half maximum (FWHM) of Mrk 279 varies with the fil- 
ter wavelength from 3".3 to 5".0, (table ||) but is always 

2 CIA is a joint development by the ESA Astrophysics Divi- 
sion and the CAM Consortium. 
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consistent with that of a point source. Monochromatic in- 
tensities were obtained by integration of the source flux 
within a circle of radius 6 pixels (18") and subtraction of 
a normalised background measured in a concentric circu- 
lar ring. The intensity of Mrk 279 in the different filters is 
listed in Table 2. The accuracy of these measurements is 
±10%. 

The ISOPHOT-S data were red uced w ith the PHOT 
Interactive Analysis (PIA; ^| Gabriel 1998] ) software pack- 
age. However, because ISOPHOT-S was operating close 
to its sensitivity limit, special reduction and calibration 
procedures had to be applied. After a change of illumina- 
tion, the responsivity of the Si:Ga photoconductors imme- 
diately jumps to an intermediate level. This initial jump is 
followed by a characteristic slow transient to the final level. 
At the faint flux limit, this time constant is extremely 
long, and in practice only the initial step is observed in 
chopped-mode. The spectral response function for this 
particular mode and flux-level was derived directly from 
observations of a faint standard star HD 132142 whose flux 
ranges from 0.15 to 2.54 Jy. The calibration star observa- 
tion was performed with the same chopper frequency and 
readout-timing as the AGN observations. The S/N of the 
ISOPHOT-S spectra was considerably enhanced by two 
additional measures: i) the 32-s integration ramps were di- 
vided into sub-ramps of 2 sec and no de-glitching (removal 
of cosmic ray hits) was performed at ramp-level it) after 
slope-fitting and de-glitching at slope-level, the maximum 
of the distribution of the slopes was determined by fitting 
a gaussian to the histogram. The resulting ISOPHOT-S 
fluxes are accurate to within ±10 %. 

The 1996 February 5 PHT-S spectrum is shown in 
Fig. 0, together with the monochromatic intensities mea- 
sured with ISOCAM on the same day. Each PHT-S data- 
point is shown with its formal error as propagated by the 
PIA software. The CAM and PHT-S fluxes agree to better 
than ±10%, providing further confidence in the reliability 
of the flux calibration. 



3. MIR flux Measurements and uncertainties 

3.1. PHT-S flux measurements and reproducibility 

For each epoch of observation the continuum flux was mea- 
sured over two different intervals, at short wavelengths 
(SW: 2.5-4.7^m) and long wavelengths (LW: 5.8-9.9 /mh). 
Table [l] lists the mean intensities over these intervals and 
their uncertainties, while the light curves are shown in 

Fig. |. 

An accurate determination of the flux uncertainties 
is essential when discussing source variability. We have 
therefore investigated the different source of errors which 
could potentially affect our PHT-S measurement. 




2xl0 13 



5xl0 13 
Rest frequency (Hz) 



10 1 



3 PIA is a joint development by the ESA Astrophysics Divi- 
sion and the ISOPHOT Consortium led by the ISOPHOT PI, 
D. Lemke, MPIA, Heidelberg. 



Fig. 1. The ISOPHOT-S spectrum from 1996 February 5 
together with the ISOCAM photometric fluxes from the 
same date. The x-axis shows the rest frequency at the 
bottom and the rest wavelength in micron at the top, both 
on a logarithmic scale 

In staring mode, the overall responsivity of PHT-S is 
known to remain stable within ± 10% (Schulz 1999). We 
have assessed the stability of the PHT-S responsivity more 
specifically at the time of each of the Mrk 279 observations 
and verified that no other systematic effects were present. 
For this purpose, two different types of calibration mea- 
surements were used as diagnostic: 

1. The detector dark current measurements, which are 
obtained immediately prior to each PHT-S observa- 
tion: throughout the campaign, the dark current sig- 
nal retained its nominal value of ~ zero V/s. We can 
therefore be confident that none of the Mrk 279 obser- 
vations suffered from detector remanence induced by 
a prior exposure to a bright source. 

2. Measurements with an internal calibration source 
which are carried-out systematically at the beginning 
of each ISO revolution: averaged over all pixels, the 
variations of detector responsivity from epoch to epoch 
are ~ ± 2%, with upper limits of 3% and 5% in the 
SW and LW range, respectively. 

A conservative upper limit of ±10% was thus adopted 
for the systematic uncertainty on the PHT-S fluxes of 
Mrk 279. Internal measurement errors were added in 
quadrature to this systematic uncertainty. The internal 
errors were computed as the dispersion about the mean 
flux in the SW and LW integration intervals, after nor- 
malization of the spectra. The purpose of the normal- 
ization is to remove the spectral curvature. Each spec- 
trum was first divided by its best-fit power-law continuum 
(F v = 8.44 lO" 14 ^- - 8 ergs^cm^Hz- 1 ; see §§) and the 
rms deviations were computed. The results show that the 
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internal errors associated with each flux measurement are 
of 2.8% and 1.7% for the 2.52-4.70 /im and 5.76-9.89 /jm 
bands, respectively. 

As a consistency check, errors were also computed by 
comparing PHT-S fluxes obtained within 30 days from 
each others. This gives a conservative error estimate since 
it assumes that there are no flux variations on time scales 
shorter than 30 days. Taking every pair of fluxes within 
30 days and measuring the error on their means, we get 
mean relative errors of 3.5% and 1.6% for the SW and 
LW bands, respectively. 



Table 1. PHT-S observation log and fluxes 

UT MJD F(2.5-4.7/t) F(5.8-9.9/t) 







(-2,450,000) 


(mJy) 


(mJy) 




(1) 


(2) 


(3) 




(4) 


1996 


Feb 5 


119 


73.5 


± 


7.6 


132 ± 13 


1996 


Mar 3 


146 


75.3 


± 


7.8 


123 ± 12 


1996 


Mar 12 


155 


65.2 


± 


6.8 


125 ± 13 


1996 


Apr 2 


176 


69.5 


± 


7.2 


129 ± 13 


1996 


Apr 27 


201 


72.7 


± 


7.6 


127 ± 13 


1996 


May 11 


215 


72.8 


± 


7.6 


128 ± 13 


1996 


May 29 


233 


73.2 


± 


7.6 


131 ± 13 


1996 


Jul 29 


294 


63.6 


± 


6.6 


127 ± 13 


1996 


Aug 12 


308 


61.8 


± 


6.4 


122 ± 12 


1996 


Aug 27 


323 


62.0 


± 


6.4 


123 ± 12 


1996 


Sep 15 


342 


60.0 


± 


6.2 


124 ± 13 


1996 


Oct 17 


374 


73.8 


± 


7.7 


128 ± 13 


1996 


Nov 1 


389 


62.6 


± 


6.5 


120 ± 12 


1996 


Nov 18 


406 


73.7 


± 


7.7 


128 ± 13 


1996 


Dec 5 


423 


70.1 


± 


7.3 


124 ± 13 


1997 


Feb 13 


493 


81.8 


± 


8.5 


128 ± 13 



Table 2. ISOCAM narrow band filter intensities 



Filter 




Range 


Flux 


FWHM 




(H 


(H 


(mJy) 


(") 


(1) 


(2) 


(3) 


(4) 


(5) 


SW1 


3.57 


3.05-4.10 


68 


3.9 


SW5 


4.25 


3.00-5.5 


60 


3.3 


LW4 


6.00 


5.50-6.50 


106 


4.1 


LW2 


6.75 


5.00-8.50 


115 


4.8 


LW5 


6.75 


6.50-7.00 


106 


4.2 


LW6 


7.75 


7.00-8.50 


120 


5.0 


LW7 


9.62 


8.50-10.7 


159 


3.8 


LW8 


11.4 


10.7-12.0 


212 


4.5 


LW3 


15.0 


12.0-18 


209 


5.0 



3.2. Comparison with ground based measurements and 
estimation of the host galaxy contribution 

Spinoglio et al. (1985) measured L-band (~ 3.5 ^m) fluxes 
of 100 ± 21mJy, 112 ± 27mJy, and 68 ± 15mJy through 
apertures of 12", 12", and 17", respectively, consistent 
with our results to within the measurement uncertainties. 

Given the spectrograph aperture (24" x 24"), the host 
galaxy of Mrk 279 could, in principle, contribute to the 
PHT-S flux. Indeed, a faint extended nebulosity is appar- 
ent in the X-band (~ 2.2 /im) image of McLeod & Rieke 
(1995). This extended flux arises from the integrated emis- 
sion of giants and supergiants in the galactic disk whose 
energy distribution is maximum at ~ 2 /xm and falls-off 
abruptly at longer wavelengths. In practice, stellar emis- 
sion will therefore make a negligible contribution to the 
PHT-S flux. Nevertheless, this was positively verified by 
comparison with ground-based data as follows: 

1. In the if-band, McLeod & Rieke (1995) estimated that 
the AGN contributes 90% of the flux within the central 
1"5 (FWHM) and 55% (35 mJy) of the total i^-band 
flux (68 mJy) integrated over the whole galaxy (i.e., 
out to a radius of 35"). Assuming "normal" near-IR 
colors (K-L = 0.22 ±0.02 mag; Clavel et al. 1989) for 
the stellar population, we derive a total host-galaxy 
flux of 18mJy at 3.5 /im. 

2. Using the -B-and i?-band nucleus-galaxy decomposi- 
tion of Granato et al. (1993), and assuming normal 
optical-to-IR (V - K = 3.22 mag; Clavel et al. 1989) 
and K — L colors, we estimate a total galaxy flux of 
20mJy and 17mJy respectively at 3.5 /im. 

These consistent estimates can be used to infer the amount 
of stellar light which enters the 24" x 24" PHT-S aperture, 
5 ± 2 mJy at 3.5 /im flux. Such a small contamination is 
within the measurement errors and can be neglected. 

The galaxy contributes 11 ± 4mJy, 14 ± 5mJy and 
11 ± 4mJy to the 12" photometric measurements in the 
J, H, and K bands respectively (Granato et al. 1993). 
These values are used in §^ to infer the intrinsic spectral 
energy distribution of the active nucleus in Mrk 279 

4. Optical observations and data reduction 

Ground-based optical observations were made in support 
of the ISO observations. Spectroscopic monitoring was 
carried out with the 1.8-m Perkins Telescope of the Ohio 
State and Ohio Wesleyan University at the Lowell Ob- 
servatory, the 1.0- m telescope of the Wise Observatory, 
the 2.6-m Shajn Telescope of the Crimean Observatory, 
and the 1.5-m Mt. Hopkins Telescope of the Harvard- 
Smithsonian Center for Astrophysics (CfA). A detailed 
log of the spectroscopic observations can be found at the 
International AGN Watch website. 

In addition CCD photometry was made on the 1.0-m 
telescope of the Wise Observatory. The flux was measured 
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using a fixed aperture of radius 7", under seeing conditions 
of 2-3". The B, V, R, and I instrumental magnitudes are 
listed in columns (3) - (6) of table [I| They have not been 
scaled to any standard system. 

The spectroscopic data were processed by the indi- 
vidual observers in standard fashion for CCD frames. 
However, the standard astronomical flux calibration tech- 
niques based on determining the instrument response 
function from observations of standard stars are only ac- 
curate for AGN spectrophotometry at about the 10% level 
even under ideal observing conditions. We used the stan- 
dard stars only for relative calibration, and employed dif- 
ferent calibration techniques for absolute calibration: The 
data from Wise Observatory were calibrated following the 
method described by Maoz et al. (1994). The data from 
Ohio, Crimean, and CFA were calibrated through scal- 
ing through [O in] A5007 flux that was measured on five 
photometric nights, at F([0m] A5007) = (1.512 ± 0.096) x 
10 -13 ergs s _1 cm~ 2 . 

On these calibrated spectra we measured the contin- 
uum flux by averaging the flux in a 10 A wide band cen- 
tered at 5105 A in the rest frame of Mrk 279 (F A (5100 A)). 
The H/3 line flux has been measured by linear interpolation 
between rest-frame wavelengths ~ 4765 A and ~ 5105 A, 
and line integration between 4770 A and 4935 A. The long- 
wavelength cutoff of this integration band misses some of 
the H/3 flux underneath [Oin] A4959, but avoids the need 
to estimate the Fen contribution to this feature and still 
gives a good representation of the H/3 variability. We did 
not correct for the narrow-line, which is expected to be 
constant. 

As the measurements made from the spectra are sub- 
ject to systematic differences between the four data sets 
used, we applied the prescriptions of Peterson et al.(1999) 
to intercalibrate the data sets, and correct for aperture 
effects. The final continuum F A (5100 A) and H/3 emission- 
line fluxes are given in table |5| The spectrophotometric 
and photometric light curves are shown in Fig. ^. Using 
the results of Granato et al. (1993), we estimate that stars 
contribute for 1.2±0.6 mJy or 20 ± 10% to the mean re- 
denning corrected 5100 A flux. 

5. The mid-IR spectrum of Mrk 279 

Figure || shows the weighted mean 2.5 - 11.7/xm spectrum 
of Mrk 279 obtained by averaging all 16 PHT-S spectra. 

The spectrum of Mrk 279 is quite similar to the mean 
Seyfert 1 spectrum obtained by Clavel et al. (2000) from 
their sample of 28 type 1 AGNs. It shows a strong con- 
tinuum, with a flux density per frequency unit that drops 
sharply with increasing frequency and only weak broad 
emission features. The continuum is well approximated by 
a power law (F v oc v a ) of spectral index a = —0.80 ± 0.05 
(Fig. ||), close to the average Seyfert 1 mid-IR index 
a = -0.84 ± 0.24 (Clavel et al. (2000)). Its flux at a fidu- 
cial wavelength of 7 /x is 103 mJy. While the broad emis- 
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Fig. 2. IR and optical light curves of Mrk 279. The top 
two panels show the IR light curves from the SW and LW 
detectors (2.52-4.70 /xm and 5.8-9.9 /xm, respectively). 
The third panel shows U-band photometric measurements 
from Wise Observatory. The forth panel shows the 5100 A 
continuum flux measured from the spectra, in units of 
10 -15 ergs s _1 enr 2 A -1 . The H/3 flux is shown in the bot- 
tom panel, in units of 10 -13 ergs s _1 cm~ 2 . The optical 
photometry, optical continuum and H/3 fluxes are avail- 
able in electronic form at the CDS. 



sion features of Polycyclic Aromatic Hydrocarbon (PAH) 
bands (Puget et al. 1985) at 3.3 /jm, 6.2 /xm, 7.7 /xm, and 
8.6 /xm are ubiquitous in many different galactic and ex- 
tragalactic line of sight, only the strongest band at 7.7 /xm 
is unambiguously detected in Mrk 279, with an intensity 
of 0.76 ± 0.12 mJy. Clavel et al. (2000) showed that PAH 
emission in AGNs originates in the interstellar medium 
(ISM) of the galaxy, whereas the mid-IR power-law con- 
tinuum arises from near nuclear dust emission in the torus. 
Because Mrk 279 is a luminous AGN, almost a quasar, the 
apparent weakness of its PAH emission can be understood 
as a contrast effect whereby a faint ISM is observed against 
a bright nucleus. The 9.7 /xm silicate absorption feature, 
conspicuous in the mid-IR spectra of starburst galaxies 
(Moorwood et al. 1996; Rigopoulou et al. 1996; Acosta- 
Pulido et al. 1996), is absent from the Mrk 279 spectrum. 



6 



M. Santos Lleo et al.: Monitoring of the mid-infrared flux of Mrk 279 with ISO 



10 



300 



,§00 



100 




-10 

CD 
I 

S 



5xl0 13 

Rest frequency (Hz) 



10' 



(http: / /tartarus. gsfc.nasa.gov/) 



-10.5 



to 

a 



Fig. 3. The weighted averaged mid-IR spectrum of Mrk 
279 (dark line), the rms deviations about the mean (light 
line) and the error (dotted line). The best-fit power law 
with index a = —0.8 is shown as a dashed line. The ex- 
pected positions of PAH and Bra features are indicated. 
The x-axis shows the rest frequency at the bottom and the 
rest wavelength, in microns, at the top, both represented 
on a logarithmic scale 



6. The spectral energy distribution 

Figure ^ shows the continuum spectral energy distribu- 
tion (SED) of Mrk 279 from the far-IR to the X-rays. The 
IRAS data points are the average of 6 pointed observations 
reported by Edelson & Malkan (1987). The near-IR data 
are from Spinoglio et al. (1985) after subtraction of the 
stellar light (§3.2). The R and B band fluxes (Granato et 
al. 1993) and the mean 5100A-flux (this paper) have also 
been corrected for the underlying galaxy contribution and 
de-redenned. The 1500 A data point represents the aver- 
age of 26 observations made with IUE (Rodriguez et al. 
1998) between 1978 and 1991. It has been corrected for 
foreground redenning using Nn = 1.6 x 10 20 cm -2 (Elvis 
et al. 1989). The large error bar reflects the strong variabil- 
ity of Mrk 279 at UV wavelengths. The EXOSAT X-ray 
data are from 1983 and 1984 and are best described in 
terms of a broken power-law (Ghosh & Soundararajape- 
rumal 1992), while the 1994 data from ASCA are mod- 
eled with a unique power law in the Tartarus Database 



Since this SED is constructed from data collected over 
~ 19 years, we caution that it may be distorted by variabil- 
ity. Variability is important in shaping the X-ray, ultra- 
violet, and optical spectrum, but is much less significant 
at longer wavelengths (see §[?]). Bearing these limitations 
in mind, it is still possible to draw some general conclu- 
sions which are not affected by flux variations at short 
wavelengths. 



■ 1 i 1 

far IR 
IRAS 



nearlR opt-UV 
ISO IUE 



I 1 1 1 1 I 1 1 

X— rays 

EXOSAT 




12 13 14 15 16 17 IB 

Log V (Hz) 

Fig. 4. The spectral energy distribution (SED) of 
Mrk 279, from the far-IR to the X rays. The IRAS data 
are represented as open squares whereas near-IR ground- 
based data are shown as stars. The PHT-S spectrum is 
shown as a dotted line, while ISOCAM fluxes are plot- 
ted as crosses where the horizontal bars indicate the filter 
range. The two open circles in the optical are the nuclear 
R- and B- band fluxes. The filled circle is the mean 5100A 
flux. Optical and near IR data have been corrected for stel- 
lar light while UV and optical data have been corrected 
for galactic redenning. The best fit power laws mid-IR and 
UV (a = —0.7) continua are also displayed. The error bar 
on the UV flux represents the rms fluctuation about the 
mean 1500A due to variability. This SED was assembled 
from data gathered over a time span of ~19 years. 



The Mrk 279 SED displays three broad maxima or 
"bumps" . The first maximum occurs in the far-IR at wave- 
lengths > 25 fim. Given the large IRAS (<~ 1') aper- 
tures and the cold color temperature of the far-IR bump, 
the 100 jum and 60 (im fluxes are probably dominated by 
cold dust from the host galaxy's ISM. Hence, the far-IR 
"bump" is most likely not related to the AGN itself. The 
second maximum is the usual "big blue bump" which dom- 
inates the SED of type 1 AGNs from the optical to the 
soft X-rays. It is usually identified as thermal emission 
from an accretion disk. In between these two maxima lies 
a third and smaller bump which extends from ~ 1 fjxa to 
~ 15-20 /urn. We tentatively identify this mid-IR bump 
as thermal emission from dust in the putative molecular 
torus and/or from dust in the NLR, as discussed below. 



7. The IR variability 

The total duration of our ISO campaign was 374 days, 
with a mean sampling interval of 24.9 days. In addition 
to the weighted-mean spectrum, (F v ), Fig. |^ displays the 
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rms spectrum (rms„), and the mean error spectrum e„. 
The latter was evaluated as 

^ = N EL (1) 

where A„ t i is the uncertainty associated with the flux Fi 
at epoch i and N = 16 is the number of epochs. We 
can exclude significant (i.e. 3cr) flux variability since the 
largest value of (rms„)/e„ is only rs 1.4. A chi-square test 
was also applied to the average PHT-S fluxes of Table |l|. 
The reduced chi-squares are < 1 which confirms that no 
significant variations of the mid-IR flux took place during 
our ~ 1-year observing campaign. 

8. Optical variability and the H/3 emission-line lag 

The relatively dense sampling of the optical light curves 
between January and July 1996 allows us to measure the 
time-delayed response, or "lag", of the H/3 emission line 
to continuum variations by cross-correlation of the light 
curves shown in Fig. 0. We used both the interpolation 
method of Gaskell & Sparke (1986) and the discrete- 
correlation function (DCF) method of Edelson & Kro- 
lik (1988), in both cases employing the specific imple- 
mentation described by White & Peterson (1994). The 
centroid of the cross-correlation is at 16.7+|'g days. Un- 
certainties were estimated using the model-independent 
FR/RSS Monte-Carlo method described by Peterson et al. 
(1998). The cross-correlation results are shown in Table || 
and Fig. |[ For comparison, Table || also lists the results 
of previous Mrk 279 monitoring campaigns by Maoz et al. 
(1990) and Stirpe et al. (1994). 



9. Discussion and concluding remarks 

The mid-IR spectrum of Mrk 279 shows a strong power- 
law continuum of spectral index —0.80 ± 0.05, with weak 
PAH emission bands and no detectable silicate ~ 9.7 /jm 
feature. The mid-IR bump of Mrk 279 extends^] from 
roughly 1.25 ^m to 15-20 /Ltm and is wider than a single 
blackbody. It peaks near 3 /im . In Fairall 9, the mid- 
IR bump most likely originates from the re-processing of 
UV and optical photons by nuclear dust (§1). We can es- 
timate the distance to the central source of the innermost 
and hottest dust grains in Mrk 279, n n , by scaling directly 
from Fairall 9 (Clavel et al. 1989). Mrk 279 is approxi- 
mately eight times less luminous than Fairall 9. Since the 
inner radius of the dust distribution is presumably con- 
trolled by sublimation, rj n should scale approximately as 
L 1 / 2 , so r; n should be a factor of ~ 2.8 smaller in Mrk 279 
than in Fairall 9, i.e. r- m w 140 ± 36 light-days. 

4 for the extension of the bump we consider the wavelength 
range AA over which uF v exceeds one third of its peak value. 




-40 -20 20 

Lag (days) 



Fig. 5. Cross-correlation function of the H/3 line inten- 
sity versus the 5100 A optical continuum flux between JD 
2,450,095.1 and JD 2,450,289.7 (Fig. g). The interpolation 
CCF is shown as a solid line and the points with error bars 
are computed using the DCF method (with a bin size of 
8 days). The H/3 line lags behind the continuum by about 
16 days. 



During the ISO campaign, the mid-IR flux did not ex- 
perience variations of amplitude larger than 10 %, the de- 
tection limit of the PHT-S instrument. Optical data con- 
temporaneous to the IR observations revealed significant 
fluctuations of the 5100 A flux with a relative rms ampli- 
tude of 9% and a ratio of the maximum to the minimum 
fluxes, R max ~ 1.64±0. 11. Any upper limit to the mid-IR 
variability in Mrk 279 has to be examined in the light of 
the UV and optical continuum variations over the same 
period of time. As noted earlier, in the dust reprocess- 
ing scenario the amplitude of the MIR flux variations will 
be reduced compared to that of the primary UV-optical 
source because of the finite propagation time of the pho- 
tons. Imagine a short (duration < 1 day) pulse of the UV- 
optical source illuminating a thin dust annulus, inclined 
by i = 10° with respect to the line of sight. The annulus 
IR response will be delayed by St = (1 — sini)n n /c and 
will last for 2(sini)ri n /c. For numerical values appropriate 
to Mrk 279, the duration of the IR reverberated pulse will 
be 38 days and its peak amplitude thereby reduced by a 
factor of order 38. Given the relatively low amplitude of 
the optical flux variations (Fig. ||), the absence of mea- 
surable variations of the MIR flux is consistent with the 
above scenario. 

Though a detailed quantitative fit with a particular 
model is beyond the scope of this paper, it is neverthe- 
less illustrative to perform a qualitative comparison of 
our data with the theoretical predictions from the torus 
model by Pier & Krolik (1992). This model predicts a mid- 
IR "bump" that is approximately 0.7 to 1 decade wide in 



8 



M. Santos Lleo et al.: Monitoring of the mid-infrared flux of Mrk 279 with ISO 



Table 3. Cross-Correlation Results 





This Work (H/3) 


Wise Obs.(H/3) 


Wise Obs.(Ha) 


LAG (Ho 


Parameter 


(1996) 


(1988) 


(1988) 


(1990) 


(1) 


(2) 


(3) 


(4) 


(5) 


Cross-correlation centroid r CGnt (days) 


16.7±S;1 


n r+25.5 
z -°-5.4 


11.6±5i 8 T 


6.8+_Z* 


Cross-correlation peak r pca k (days) 


2l + _l 


q + 28 
J -5 


11+n 
ii-n 


0+29 
A -3 


Peak correlation coefficient r max 


0.769 


0.799 


0.793 


0.795 


Continuum rms fractional variability F val 


0.093 


0.071 


0.071 


0.126 


Mean time between observations (days) 


2.3 


4.1 


4.1 


6.0 


Duration of campaign (days) 


195 


156 


156 


152 



wavelengths, in agreement with the Mrk 279 observations. 
In the Pier & Krolik (1992) model, the torus emission is 
expected to peak at a wavelength A poa k that depends pri- 
marily on the flux illuminating the torus inner surface and 
its inclination angle i with respect to the line of sight. The 
relatively high color temperature implied by A poa k ~ 3 /urn 
constrains the inclination to be small (cosi > 0.75). The 
absence of silicate absorption also rules-out very optically 
thick models and constrains the vertical column density 
at the torus inner edge, Nh < 10 24 cm~ 2 . Comparison of 
Fig. § with Fig. 5 of Pier & Krolik (1992) also suggests a 
moderately thick torus, with r- m /h = 0.3. 

The delay AT of H/3 w.r.t. the optical continuum was 
16.7l 5 ; 6 days during this campaign. Comparison with the 
results from previous monitoring campaigns (see Table [}]) 
does not reveal any significant change of AT over a time 
span of ~8 years. In other words, we find no evidence for 
a secular change in the structure of the BLR in Mrk 279. 
Equating c x AT with the emissivity weighted radius 
Rblr of the H/3 emitting region, one sees that r- m is about 
8 times larger than Rblr- In other words, the BLR lies 
well within the dust evaporation radius. 
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Table 4. Optical Photometry from Wise Observatory 



UT 


Julian Date 




B* 


V* 


R* 


I* 






Date 


(-2,450,000) 




(mag) 


(mag) 


(mag) 


(mag 









(1) 


(2) 




(3) 


(4) 


(5) 


(6) 






1995 


Dec 28 


80.6 


13. 


.838 ±0.020 


12.533 ±0.016 


12.461 ±0.012 










1996 


Jan 29 


112.4 


14. 


111 ±0.040 


12.702 ±0.018 


12.578 ±0.010 


12.771 


± 





,022 


1996 


Mar 27 


170.6 


14. 


.097 ±0.015 


12.706 ± 0.022 


12.588 ±0.016 










1996 


Mar 28 


171.5 


14. 


.070 ±0.021 


12.715 ±0.022 


12.605 ± 0.012 










1996 


Mar 29 


172.5 


14. 


151 ±0.013 


12.805 ± 0.005 


12.699 ± 0.006 










1996 


Mar 31 


174.5 


14. 


110 ±0.024 


12.738 ±0.021 


12.631 ± 0.023 


12.829 


± 


0. 


.023 


1996 


May 2 


206.5 


14. 


517 ±0.015 


13.243 ±0.014 


13.130 ±0.012 


13.426 


± 





026 


1996 


Jul 1 


266.3 


14. 


377 ± 0.042 


12.697 ±0.056 


12.733 ± 0.098 










1996 


Jul 2 


267.3 


14, 


417 ±0.030 




12.762 ±0.051 


12.964 


± 


0, 


.080 


1996 


Aug 6 


302.2 


14, 


111 ±0.017 


12.745 ±0.014 


12.626 ±0.017 


12.831 


± 


0. 


.018 


1996 


Dec 19 


437.6 


14, 


245 ± 0.045 


12.810 ±0.037 


12.659 ± 0.022 


12.837 


± 


0. 


.012 


1997 


Feb 18 


498.6 


14, 


192 ±0.013 


12.787 ±0.008 


12.673 ±0.019 


12.856 


± 





.024 


1997 


Mar 19 


527.5 


14, 


179 ±0.019 


12.769 ±0.014 


12.661 ± 0.009 


12.842 


± 


0. 


.021 


1997 


May 10 


579.4 


14, 


090 ± 0.027 


12.684 ± 0.014 


12.576 ± 0.009 


12.776 


± 


0. 


.013 


1997 


Jun 9 


609.4 


14, 


088 ± 0.021 


12.701 ±0.013 


12.604 ±0.012 










1997 


Jul 25 


655.3 


14. 


005 ± 0.028 


12.652 ±0.016 


12.550 ± 0.021 


12.770 


± 


0. 


.013 


1997 


Aug 30 


691.3 


14. 


118 ±0.030 


12.705 ±0.010 


12.578 ± 0.020 











★ Instrumental magnitudes, not scaled to any standard system 
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Table 5. Optical Continuum and H/3 fluxes 



Julian Date 


F A (5100 A) 


F(R/3) 


(-2,450,000) 


(a) 


(b) 


(1) 


(2) 


(3) 


95.1 


9.39 ±0.40 


6.56 ±0.28 


96.0 


9.67 ±0.41 


6.57 ±0.28 


97.0 


8.59 ±0.37 


6.10 ±0.26 


98.1 


9.09 ±0.39 


6.30 ±0.27 


104.0 


7.96 ±0.35 


6.36 ±0.28 


108.0 


7.84 ± 0.35 


6.32 ±0.27 


109.0 


7.37 ±0.33 


6.18 ±0.27 


110.0 


7.74 ± 0.35 


6.25 ±0.27 


110.4 


7.29 ±0.22 


6.58 ±0.20 


111.0 


7.57 ±0.34 


6.38 ±0.28 


112.1 


7.54 ± 0.34 


6.07 ±0.26 


119.0 


7.20 ±0.14 


6.66 ±0.13 


123.0 


6.99 ±0.14 


6.26 ±0.12 


124.0 


7.08 ±0.32 


6.19 ±0.27 


128.0 


7.08 ±0.32 


5.54 ±0.24 


129.0 


7.46 ± 0.34 


5.86 ±0.25 


130.0 


7.12 ±0.33 


5.63 ±0.24 


131.0 


7.53 ±0.34 


5.90 ±0.26 


132.0 


7.51 ±0.34 


5.69 ±0.25 


134.0 


7.18 ±0.33 


5.55 ±0.24 


134.9 


7.25 ±0.33 


5.75 ±0.25 


136.9 


7.42 ±0.15 


6.07 ±0.12 


143.0 


6.99 ±0.14 


5.62 ±0.11 


149.9 


6.80 ±0.14 


6.05 ±0.12 


156.5 


7.11 ±0.21 


6.02 ±0.18 


164.9 


7.11 ±0.14 


5.90 ±0.12 


168.9 


7.34 ±0.14 


5.83 ±0.11 


170.0 


7.42 ± 0.34 


5.58 ±0.24 


171.0 


7.51 ±0.34 


5.88 ±0.25 


175.9 


7.27 ±0.14 


6.07 ±0.12 


185.9 


6.73 ±0.12 


6.24 ±0.11 


188.0 


6.51 ±0.31 


5.79 ±0.25 


189.0 


6.07 ±0.29 


5.47 ±0.24 


190.0 


6.46 ±0.31 


5.94 ±0.26 


191.0 


6.61 ±0.31 


6.12 ±0.26 


191.9 


6.31 ±0.12 


5.96 ±0.11 


193.0 


6.62 ±0.31 


5.95 ±0.26 


194.0 


6.11 ±0.29 


6.09 ±0.26 


198.8 


6.07 ±0.11 


5.82 ±0.11 


199.8 


6.53 ±0.31 


5.97 ±0.26 


206.8 


6.48 ±0.13 


5.47 ±0.11 


211.8 


6.55 ±0.31 


5.70 ±0.25 


212.9 


6.41 ±0.12 


5.47 ±0.10 


213.9 


6.79 ±0.32 


5.60 ±0.24 


215.4 


6.40 ±0.19 


5.26 ±0.16 


218.6 


7.16 ±0.33 


5.24 ±0.23 


219.6 


7.11 ±0.33 


5.77 ±0.25 


220.9 


7.06 ±0.14 


5.33 ±0.11 


221.6 


6.95 ±0.32 


5.09 ±0.22 


222.7 


6.95 ±0.32 


5.26 ±0.23 


223.6 


7.26 ±0.33 


5.49 ±0.24 


224.6 


6.95 ± 0.32 


5.44 ± 0.24 



(a) 10~ 15 ergs s" 1 cm" 2 A -1 

(b) 10~ 13 ergs s" 1 cm" 2 



Table 5. Optical Continuum and H/3 fluxes (cont.) 



Julian Date 


F A (5100 A) 


F(H/3) 


(-2,450,000) 


(a) 


(b) 


(1) 


(2) 


(3) 


225.6 


6.88 ± 0.32 


5.33 ±0.23 


226.6 


6.87 ±0.32 


5.24 ±0.23 


227.8 


6.84 ±0.14 


5.43 ±0.11 


231.6 


7.06 ± 0.32 


5.69 ±0.25 


233.8 


6.82 ±0.14 


5.57 ±0.11 


238.3 


6.51 ±0.19 


5.82 ±0.17 


240.6 


6.70 ± 0.31 


6.04 ± 0.26 


241.8 


6.67 ±0.13 


5.69 ±0.11 


242.6 


6.82 ± 0.32 


6.20 ±0.27 


248.8 


6.71 ±0.13 


5.47 ±0.11 


249.7 


6.30 ± 0.30 


5.84 ±0.25 


250.7 


6.23 ± 0.30 


5.44 ± 0.24 


251.7 


6.84 ± 0.32 


5.92 ±0.26 


253.6 


5.85 ± 0.28 


5.63 ±0.24 


253.8 


6.01 ±0.12 


5.55 ±0.11 


254.6 


5.90 ± 0.29 


5.84 ±0.25 


256.6 


6.43 ± 0.30 


5.64 ±0.25 


257.7 


6.63 ± 0.31 


5.80 ±0.25 


258.6 


6.24 ± 0.30 


5.30 ±0.23 


259.6 


6.04 ± 0.29 


5.25 ±0.23 


262.8 


6.44 ±0.13 


5.35 ±0.11 


270.7 


6.20 ± 0.30 


5.27 ±0.23 


272.3 


6.83 ± 0.20 


5.07 ±0.15 


280.6 


6.62 ±0.31 


5.13 ±0.22 


281.7 


6.80 ± 0.32 


5.53 ±0.24 


282.7 


6.05 ± 0.29 


5.01 ±0.22 


283.7 


6.07 ±0.29 


5.19 ±0.22 


284.7 


6.51 ±0.31 


5.36 ±0.23 


286.7 


6.86 ± 0.32 


5.28 ±0.23 


287.6 


6.54 ±0.31 


4.96 ± 0.22 


288.7 


6.94 ± 0.32 


5.30 ±0.23 


289.7 


6.46 ± 0.31 


5.04 ±0.22 


460.1 


6.46 ±0.31 


5.54 ±0.24 


465.0 


6.44 ± 0.30 


5.55 ±0.24 


480.1 


5.98 ± 0.29 


5.66 ±0.25 


481.1 


6.39 ± 0.30 


5.46 ± 0.24 


482.1 


6.36 ± 0.30 


5.75 ±0.25 


485.9 


6.51 ±0.31 


5.84 ±0.25 


487.0 


6.43 ± 0.30 


5.94 ±0.26 


489.0 


6.58 ± 0.31 


5.79 ±0.25 


489.9 


6.79 ± 0.32 


5.96 ±0.26 


491.0 


6.55 ±0.31 


5.49 ± 0.24 


494.0 


6.30 ± 0.30 


5.42 ± 0.23 


510.0 


6.66 ±0.31 


5.60 ±0.24 


519.8 


6.61 ±0.31 


5.30 ±0.23 


546.4 


7.42 ± 0.22 


4.88 ±0.15 


547.8 


7.61 ±0.34 


5.14 ±0.22 


548.8 


7.43 ± 0.34 


5.21 ±0.23 


549.8 


7.64 ± 0.34 


5.39 ±0.23 


550.8 


7.64 ± 0.34 


5.59 ±0.24 


551.8 


7.15 ±0.33 


5.24 ±0.23 


552.8 


7.49 ± 0.34 


5.55 ±0.24 


569.8 


7.55 ± 0.34 


5.43 ±0.23 



(a) 10 -16 ergs s" 1 cm" 2 A" 1 

(b) 10~ 13 ergs s" 1 cm" 2 



12 



M. Santos Lleo et al.: Monitoring of the mid-infrared flux of Mrk 279 with ISO 



Table 5. Optical Continuum and H/3 fluxes (cont.) 



Julian Date 


F x (5100 A) 




F(H/3) 


(-2,450,000) 


(a) 




(b) 


(1) 


(2) 




(3) 


575.4 


7.78 ±0.23 


5 


.85 ±0.17 


578.8 


7.83 ±0.35 


5 


.37 ±0.23 


580.5 


7.39 ±0.22 


5 


.46 ±0.16 


581.8 


7.22 ±0.33 


5 


.52 ±0.24 


582.6 


6.96 ±0.32 


5 


.57 ±0.24 


597.7 


7.08 ±0.32 


5 


.56 ±0.24 


602.6 


0.36 ± 0.30 


5 


A A 1 A on 

.00 ± 0.22 


604.6 


6.68 ±0.31 


5 


.22 ±0.23 


605.7 


7.11 ±0.33 


4 


.93 ±0.21 


608.7 


7.41 ±0.34 


5 


.09 ±0.22 


609.6 


7.34 ± 0.33 


5 


.35 ±0.23 


610.6 


6.88 ±0.32 


4 


.99 ±0.22 


627.6 


7.45 ± 0.34 


5 


.35 ±0.23 


628.6 


7.58 ±0.34 


5 


.51 ±0.24 


629.7 


8.26 ±0.36 


5 


.83 ±0.25 


630.7 


7.92 ±0.35 


5 


.98 ±0.26 


(a) 10~ 15 ergs s" 


^cm^A" 1 







(b) 10~ 13 ergs s" 1 cm- 2 



